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Abstract. Brownian motions on a metric graph are defined, their Feller property 
is proved, and their generators are characterized. This yields a version of Feller's 
theorem for metric graphs. 



1. Introduction 



In his pioneering articles |[9- 11 1, Feller raised the problem of characterizing and 



constructing all Brownian motions on a finite or on a semi-infinite interval. In the 
sequel this problem stimulated very important research in the field of stochastic pro- 
cesses, and the problem of constructing all such Brownian motions found a complete 
solution 1 ,17, ,18J via the combination of the theory of the local time of Brownian 
motion pl| , and the theory of (strong) Markov processes @|4]^ 16 1. 



On the other hand, there is a growing interest in metric graphs, that is, piecewise 
linear varieties where the vertices may be viewed as singularities. Metric graphs arise 
naturally as models in many domains, such as physics, chemistry, computer science 
and engineering to mention just a few — we refer the interested reader to [30] for a 
review of such models and for further references. 

Therefore it is natural to extend Feller's problem to metric graphs. The present 
paper is the first in a series of three articles p4p5| (together with a more pedagogical 
one [23], in which the well-known classical cases of finite and semi-infinite intervals 
are revisited) on the characterization and the construction of all Brownian motions on 
metric graphs. Stochastic processes, in particular Brownian motions and diffusions, 
on locally one-dimensional structures, notably on graphs and networks, have already 
been studied in a number of articles, out of which we want to mention |]T][3j[8j[T2]-[T5j 
29J in this context. 

Heuristically, a metric graph {Q, d) can be thought of as the union of a collection 
of finite or semi-infinite closed intervals which are glued together at some of their 
endpoints which form the vertices of the graph, while the intervals are its edges. The 
metric d is then defined in the canonical way as the length of a shortest path between 
two points along the edges, and the length along the edges is measured as for usual 



intervals. For a more formal definition of metric graphs see section 2.1 We will 



only consider finite graphs, that is, those for which the sets of vertices and edges are 
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finite. For the definition of a Brownian motion on the metric graph {Q , d) we take a 
standpoint similar to the one of Knight |22| for the semi-line or a finite interval: It is 
a strong Markov process with cadlag paths which are continuous up to the lifetime, 
and which up to the first passage time at a vertex is a standard Brownian motion on 
the edge where it started. For the formal definition, cf. section [3] 

The crucial problem is then to characterize the behaviour of the stochastic process 
when it reaches one of the vertices of the graph Q, or in other words, the characteriza- 
tion of the boundary conditions at the vertices of the Laplace operator which generates 
the stochastic process. We want to mention in passing that in an L^-setting all bound- 
ary conditions for Laplace operators on Q which make them self-adjoint operators 



have been characterized in |26 27 1. In this series of papers we shall work with the 



Banach space of continuous functions on Q which vanish at infinity. The main result 



of the present paper is Feller's theorem for metric graphs (cf. theorem 5.3 1. Roughly 
speaking it states that all possible boundary conditions are local boundary conditions 
of Wentzell type, i.e., linear combinations of the value of the function with its first 
and second (directional) derivatives at each vertex, subject to certain conditions on 
the coefficients. We want to emphasize here, that the fact that we only obtain local 
boundary conditions is due to the assumption that the paths of the Brownian motion 
have no jumps during their lifetime. On a more technical level this assumption entails 
that we deal with Feller processes — as is proved in section|4] — which is of consider- 
able advantage when we prove the strong Markov property of the processes which we 



will construct in the follow-up papers |24 25 1. On the other hand, in a forthcoming 



work we shall embed the situation into the larger framework of Ray processes, and 
there we shall deal also with non-local boundary conditions, allowing the processes 
to have jumps from the vertices into Q in addition to the jumps to the cemetery point. 

The paper is organized in the following way. In section [2] we recall the pertinent 
notions of (finite) metric graphs, and of strong Markov processes on metric graphs, at 
the same time setting up our notation. In section[3]Brownian motions on metric graphs 
are defined, and some consequences of this definition are discussed. The proof of the 
statement that Brownian motions on metric graphs are Feller processes is given in 
section|4] Finally, we state and prove Feller's theorem for metric graphs in section[5] 
In appendix |A] we give a short account of Feller semigroups in a form which we find 
especially convenient for the purposes of the present paper, but which we could not 
find in this form elsewhere. 

The contents of the other two papers in this series are as follows. In the arti- 
cle p4| all Brownian motions on single vertex graphs (roughly speaking, n semi 



lines [0, +oo) glued together at the origin) are constructed, and in the article |25| 
these Brownian motions are pieced together pathwise to yield all possible Brownian 
motions on a general metric graph. 
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2. Preliminary Dehnitions and Notations 

2.1. Metric Graphs. Throughout this paper we consider a fixed finite metric graph 
{Q, d). That is, is a quadruple {V,I, £, d), where F is a finite set of vertices, X is 
a finite set of internal edges, <S is a finite set of external edges, and 3 is a map from 
the set £ = T U f of edges into (y x y) U V , which maps an internal edge i G X 
to an ordered pair (d^ (i\d^ {%)) G V oi vertices, called the initial and final 
vertex of i, while e £ £ is mapped to d{e) G V, called the initial vertex of e. Every 
edge / G £ is assumed to be isometrically isomorphic to an interval namely every 
e G is in one-to-one correspondence with the half line [0, +oo), while for every 
i £ I there exists pi > so that i is isomorphic to [0, pj]. Under these isomorphisms, 
for i G X we have that d-{i) corresponds to while d^{i) corresponds to pi, and for 
e £ £, the vertex d{e) corresponds to 0. pi is called the length of the internal edge 
i G X. Moreover, we suppose that the edges are sets with an ordering as induced by 
the isomorphisms mentioned above. C{v) = {Z G £, u G d{l)} is the set of edges 
incident with v. For notational simplicity we will also use d{l), I £ C, to denote the 
set consisting of d^{l) and d^{l) if Z G X, and of d{l) if I £ £. 

In [28] the standard notion of a walk on a graph (e.g., [19]) has been generalized 
to graphs of the above type, and therefore we obtain in a natural way a metric d 
on Q as the minimal length of all walks leading from one point to another, where the 
length is measured along the edges as induced by the isometry with the corresponding 
intervals. 

In the sequel it will be convenient — and without danger of confusion — to identify 
the abstract graph Q with its isomorphic geometric graph (e.g., [[T9)). In other words, 
we also consider ^ as a union of the intervals corresponding to the edges, subject 
to the equivalence relation defined by the combinatorial structure of the graph which 
identifies those endpoints of the intervals which correspond to vertices to which the 
respective edges are incident. Similarly, we shall often identify the edges with the 
intervals they are isomorphic to. 

Clearly, {Q, d) is a complete, separable metric space, and hence it is a Polish space. 
The Borel cr-algebra of Q is denoted by B{Q). We write Bt{C) for the open ball with 
radius r > and center ^ G ^. 

For I £ C,r denotes the open interior of /, i.e., /° is the subset of / being isomor- 
phic to (0, Pi) if / = i G X, and to (0, +oo) ifl £ £. We setQ° = Q\V to be the 
interior of Q, and hence Q° is the pairwise disjoint union of the open edges 1° , I £ C. 
In particular, every ^ G ^° is in one-to-one correspondence with its local coordinate 
(/, x),l £ C, X £ Ii, and we may and will write ^ = (/, x). 

Assume that / is a real valued function on Q. Then / is in one-to-one correspon- 
dence with the family of functions {fi, I £ C) where // is the restriction of / to the 
edge I £ C. (Of course, if v £ V is a. vertex with which the edges I, I' £ C are 
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incident, then we have to have fi{v) = fi'{v).) Sometimes it will also be convenient 
to write fi{x) instead of f{C), for £, ^ G° having local coordinates (/, x). 

The space of real valued, bounded measurable functions will be denoted by B{Q), 
while Co{Q) denotes the space of continuous real valued functions on Q which vanish 
at infinity. Both spaces are equipped with the sup-norm, denoted by || • ||. {B{Q), \\ ■ ||) 
and (Co(t/), || • ||) are Banach spaces, the latter being separable. 

2.2. Markov Processes on Metric Graphs. Let {Q, d) be a metric graph as in the 
previous subsection. Furthermore let {^,A) be a measurable space, equipped with 
a family (P^, ^ G ^) of probability measures. The expectation with respect to P^, 
i & Q, will be denoted by E^{-). We will say that a statement holds almost surely 
(a.s.), if for all ^ G ^ the statement holds almost surely with respect to P^. 

Let A be a point not in Q which we will view as a cemetery point. By we 
denote the union Q U {A}, where A is adjoined to Q as an isolated point. We define 
the (j-algebra B{Q^) on as the obvious minimal extension of B{Q). All real 
valued functions / on ^ are understood as being extended to with /(A) = 0. 

We consider a ^^-valued normal homogeneous Markov process X = {Xt, t > 0) 
on (17, A) relative to a filtration F = {Ft, i > 0) in ^ with cadlag paths. Throughout, 
we will suppose that the filtration F is right continuous and complete for the family 
(P^, ^ G Q), that is, for all t > ^, Ft = r\e>oFt+e, and Fq contains all subsets of Q. 
which are negligible for all P^, ^ ^ Q. Also A is a cemetery state for X, i.e., almost 
surely = A, s > 0, entails Xt = /S. for all t > s. The lifetime of X is defined 
by C = inf{i >0,Xt = A}. 

As in | [33| we assume that the transition probabilities of X are given in terms of a 
transition function P = {Pt, t > 0), i.e., 

P^iXt eC) = PtiC, C), e e t > 0, C G B{g^). 

In particular, for alH > 0, C G 8(0'^), the mapping (, ^ P^ {Xt G C) is measurable. 
In terms of the transition function P the Markov property of X can be written as 
follows: 

P^{Xt £C\Fs) = Pt-siC, C), < s < t, C G BiG^), e G 

It will be convenient and there is no loss of generality to assume the existence of a 
shift operator 9 : IR+ x $7 — 0, (s, cj) i— )• 9s{uj), so that a.s. for all t, s > 0, 

(2.1) Xto9s = Xt+s- 

Thus if X is a strong Markov process with respect to F, its strong Markov property 
can be expressed in the following way. Let S be an J^-stopping time, and as is usual 
denote cr-algebra of the past of S by Fs- If C S Q, and Z is a positive or bounded 
random variable, then 

(2.2) E^{Zo9s\Fs) =Exs(.Z), 

holds P^-a.s. on the set {Xs / A} = {5'<C}- 

For a subset A of C/ we shall denote its hitting time by X by Ha, 

Ha = mf{t >0,Xt£ A} 
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and if A = {^}, ^ £ Q, we simply write H^. 

We shall occasionally take the liberty to write X{t) for Xt, t G R+, or H{A) for 
Ha, a C Q, whenever it is typographically more convenient. 

The semigroup U = {Ut, t > 0) associated with X and acting linearly on B{Q) is 
defined by 

(2.3) UtfiC) = E^ifiXt)) = I f{7i)m,d7]), 

Jg 

for / G B{g), ^gG. The bound \\Utf\\ < \\f\\ holds for all / G B{g), t > 0. The 
resolvent R = {R\, A > 0) associated with X, and acting on is defined by 

/•oo 

(2.4) RxfiC) = / e-'^'UtfiO dt, A > 0, 

Jo 

and satisfies 

(2.5) ||i?A/||<^||/||. 

We shall denote the restrictions of the semigroup U and the resolvent R to Co(^) by 
the same symbols. 

Assume that X is a strong Markov process with respect to the filtration T. A direct 
consequence of the strong Markov property is the first passage time formula for the 



resolvent (e.g., ||32| or |[18|): Let ^ G ^, / G B{g), A > 0, and assume that 5 is a 
J^-stopping time which is P^-a.s. finite. Then 

(2.6) RxfiO = {Jj e-^'fiXt) dt) + {e'^^ Rxf{Xs)) 

holds true. 

We shall often have occasion to use a standard Brownian family on the real line M 
as a family of reference processes: Let {Q', A') denote another measurable space with 
a family {Qx, x G M) of probability measures, and for every x G M, {Bt, t G M_|_) 
a standard Brownian motion on M starting Q^^-a.s. in x. Expectations with respect 
to Qx will be denoted by Ex{ ■ )■ The Brownian family is equipped with a filtration 
denoted by J' = {Jt, t > 0), and throughout we assume — as we may — that J 
is right continuous and complete for the family {Qx, x G M). (For example, we can 
always consider the natural filtration generated by {Bt, t G and then choose its 



universal augmentation, e.g., 1 33 Chapter in.2] or 1 21 Chapter 2.7].) denotes 
the hitting time of the set ^4 C M by B, and as above we simply write for H^^y 

and occasionally B{t) for Bt, H^{A) for H^, AcR. 

Consider an edge I G C, and the interval Ii C M+ that / is isomorphic to. We set 
d{Ii) = {0, pi} iflGl, and d{Ii) = {0} if / G For t G M+ set 

Bi = B{tAHi^j^^) 

where s At = mm{s,t}, s, t G M+, i.e., is a Brownian motion on M with 
absorption in the endpoint(s) of Ii. Under the family {Qx, x G Ii) we call this process 
the absorbed Brownian motion on Ii. 
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3. Definition of Brownian Motions on a Metric Graph 



In analogy with p2| Chapter 6] we define a Brownian motion on a metric graph Q 
as follows. 

Definition 3.1. A Brownian motion on Q is a normal strong Markov process X = 
{Xt, t > 0) with state space and lifetime Q. The sample paths of X are right 
continuous with left limits in Q, and they are continuous on [0, C)- Furthermore, for 
every ^ G ^° with local coordinates {l,x), I G C, x G the process X"''" = 
(yX{t A Hy), t > O), with start in ^ and absorption in the set of vertices V, is 
equivalent to an absorbed Brownian motion on Ii with start in x. 



Remark 3.2. According to our convention of subsection 2.2 we consider a Brownian 
motion X on ^ as a strong Markov process with respect to a filtration F = {Ft , t > 
0) which is right continuous and complete. 

Remark 3.3. Suppose that ^ £ I, I € C, then Hy = Hq(j^^, because paths starting at ^ 
hit the set d{l) before any other vertex due to the continuity assumption. 



Consider a Brownian motion X onQ. Let us state the last condition in definition |3.1| 
more explicitly. Fix I G £, and ^ G /° with local coordinates (Z,x). Then for all 
n G N, ti, . . . , £ with ti < t2 < ■ • • < tn, and all Ai, An in the Borel 
cr-algebra B{1) oil, 

P^{Xt, e Ai,...,Xt,^e An, t„ < Hy) 
(3.1) 5 . 

= Qx[Bti G Ai, . . . ,Bt„ G An,tn < 

For simplicity we have identified the set Ai C / with its isomorphic image in //. 
Observe that in particular under P^, the stopping time Hy = i^a(z) has the same law 
as Hq,j . under Q^, and especially we get that for all G ^, P^{Hy < +oo) = 1. 



It follows from definition 3.1 that any discontinuity of the paths of X can only 
occur at the vertices of Q, and it consists in a jump to the cemetery state A. Hence 
if the process starts in G Q°, it cannot reach the cemetery state A before hitting V. 
On the other hand, if the process starts in w G F , P^-a.s. it cannot jump right away to 
A, because this would contradict the right continuity of the paths and the requirement 
p^(Xo = = 1 for all C G a. In particular, we have for all (,£G,P^{C> Hy) = 1. 

For the following discussion we assume that the process X starts at a vertex v G 
V, and consider the exit time from v, i.e., the stopping time Sy = H{Q°). It is 
well known (e.g., ||7]22j 33 1) that because of the strong Markov property of X, is 



under P^ exponentially distributed with a rate fiy G [0, +oo\. Thus there are three 
possibilities: 

Case = 0: In this case the process stays at v forever, i.e., u is a trap, and the 
process is given by {X{t A H^), t G 

Case < /3y < +oo: In this case the process stays at v P^-a.s. for a strictly positive, 
finite moment of time, i.e., v is exponentially holding. It is well known (cf., e.g., 1221 



p. 154], |33 p. 104, Prop. 3.13]) that then the process has to leave V by a jump, and by 
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our assumption of path continuity on [0, ?]), the process has to jump to the cemetery 
A. 

Case (3y = +00: In this case the X leaves the vertex v immediately, and it begins a 
Brownian excursion into one of the edges incident with the vertex v. 

4. Feller Property 

In this section we prove that the semigroup U associated with every Brownian mo- 
tion on Q has the Feller property (see, e.g., [22][33|, or definition A. 1 in appendix [A]). 
We begin with the following simple lemma. 

Lemma 4.1. Assume that ^ G Q°. Then P^—a.s. converges to zero, whenever rj 
increases or decreases to ^ along the edge to which ^ belongs. 

Proof. First we notice that because up to time C the paths of X are continuous, rj 1— )• 
is Pg-a.s. monotone decreasing as t] increases or decreases to ^. Therefore it is 
enough to show that Hrj converges to zero in P^-probability. 

Let ^ £ 1° , I ^ C, with local coordinates (Z, x), x G 1°. Fix e > small enough so 
that (/, X ± e) € 1°. Without loss of generality we may assume that t]) < e. We 
consider first the case where | ^, i.e., for the local coordinates {I, y), y G Ii, of rj 
we have y \^ x. Let 5 > 0, and write 

(4.1) P^{Hr, >S)= P^{H^ > 6,H(^i^,_,) > H^) + P^{Hr, > 5, < Hr,). 
We estimate the second probability on the right hand side from above by 

But this is the probability of the event that the process leaves the set on I which in 
local coordinates is the interval [x — e, y]) at the end point with local coordinates 
{l,x — e). Therefore this is an event which happens before the process hits a vertex, 
and hence this probability is equal to the corresponding one for a standard Brownian 
motion (e.g., y_8j Problem 6, p. 29]): 

which converges to zero as r/ | ^. Similarly, the first probability on the right hand 
side of equation ( |4. 1| ) is equal to 



00 



fi'y - X 

err 



sec- 



where we used the well-known density of Hy under Q^, e.g., |34 p. 292], 1 18 
tion 1.7], or pJ\ Proposition 2.6.19]. Clearly, the last expression converges to zero as 
y i X, i.e., asT] I 
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The case f C is treated in an analogous way. □ 
Lemma 4.2. Let \ > Q, v &V, and suppose that I G C{v). Then 

(4.2) lim EJe~^^'') = l. 

Proof. Fix e > in such a way that we have for every v' & V,v' ^ v, d{v, v') > e. 
We may assume that d{v, rf) < e. Set 

(4.3) H,,, = H{B,{vy), 

where the superscript c denotes the complement of a set. Write 

(4.4) 1 - Erj (e-^^") =E^{1- e-^^" ; < H,^,) + (l - e"^^" ; > H,,,) , 
with the notation 

(4.5) Er,{Z;C) = Er,{Zlc) 

for positive or P^-integrable random variables Z, and C ^ A. Consider the case 
where the vertex v corresponds to the point with local coordinates (/,0), the case 
where v corresponds to {I, pi) can be dealt with by an analogous argument. Let rj 
have local coordinates {l,y), < y < e. The second term on the right hand side of 



equation ( |4.4[ ) is less or equal to 

Pr){Hv > Hv,e) = Qy{HQ > Hf) 

= y 

e 

which converges to zero with y i 0, i.e., with r] ^ v. On the other hand 

E,{1- e-'^^ ; < H,,,) = E^ {l - e'^^o ; < ^i.) 

= 1 _ e-^^y, 

where we used the well-known formula for the Laplace transform of the density of 



under Qy (e.g., 1 18 p. 26, eq. 5]). Obviously this converges to zero as y | 0, i.e.. 



as — 7- w. □ 
Theorem 4.3. Every Brownian motion on Q is a Feller process. 

m 



Proof. The proof is based on the first passage time formula p.6[ ). By theorem A.3 
appendix [A| it suffices to show that for all A > 0, i?A maps Cq{Q) into itself, and that 
for all / G CQ{g), CgG, Utf{0 converges to /(O as t i 0. 

First we show that for every A > 0, i?A maps Co{Q) into itself. Assume that 



/ € Co{Q). Consider the case G Q°. Then it follows from lemma 4.1 as in 1 18 
Section 3.6] that Rxf is continuous at ^. Consider now the case ^ = v ^ V, let I 
belong to the set C{v) of edges incident with v, and let t] £ I. Note that P^-a.s. is 
finite (cf. section [3]l. Therefore we can employ equation ( |2.6[ ) with S = Hy-. 

Rxfiv) = e~^'f{Xt) dt) + £^(6-^^") Rxf{v). 
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Using ( |2.5[ ) we estimate in the following way 



Rxf{v)-Rxf{v) < 



e-^'f{Xt) dt 







< 



A 



1 



-XHy 



\Rxfiv)\ 



By lemma 4.2 this term converges to zero as r] converges to v along /. Therefore R\f 
is also continuous at v. 

Next we prove that for all A > 0, / G Co{Q), R\f vanishes at infinity. If Q has no 
external edges there is nothing to prove, and so we assume that e G <S is an external 
edge of Q which is incident with the vertex v ^ V , d{e) = {v}. Let A, e > be 
given. We choose ri > large enough so that for all ^ G e with d{v, ^) > ri we have 
1/(^)1 < eA/2. Choose r2 > r\, and consider ^ G e with d{v,^) > r2. Denote by 
^1 the point on e which has distance ri to v. Then we have that P^-a.s., i/^^ < Hy, 
and consequently P^{H^-^ < +oo) = 1. Hence we can use the first passage time 
formula ( |2.6| ) in the form 

Rxf iO = 



H, 



e-^'f{Xt) dt 



Rxfiii 



For t G we have d{v,Xt) > ri, and therefore the absolute value of the 

first term on the right hand side is bounded from above by e/2. For the second term 
we can compute the expectation as for the corresponding expression of the standard 
Brownian motion B on M, and we obtain (again with ( |2.5| )) 



Ed 



Rxf i^i 



< e~ 



2A{r2-ri) 



A 



Now choose r2 large enough so as to make the last term less than e/2, and we are 
done. Thus we have shown that for every X > 0, Rx maps Co{Q) into itself. 
Finally, consider for f G CoiG), C G, t > 0, 

Utf{0 = Edf{Xt)). 

By definition, X has right continuous sample paths, and P^{Xq = ^) = 1. Since 
/ is continuous and bounded, an application of the dominated convergence theorem 
shows that Utf{C) converges to /(^) as t decreases to 0. □ 

5. Generators AND Feller's Theorem 
Let Vc denote the subset of y x £ given by 

Vc = {{v,l), veV wAle C{v)]. 
We shall also write vi for (w , /) G Vc- We remark in passing that 

\Vc\ = \£\+2\I\. 
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Consider a real valued function / on Q, let u G y and let I G C{v) be an edge 
incident with v. We define the directional derivative of f at v in direction I G C{v) 
as follows: 



(5.1) f'ivi) 



lim /'(O) if u is an initial vertex of 
— lim f'(£), if u is a final vertex of I, 



whenever the corresponding limit on the right hand side exists. Geometrically this 
directional derivative is just the inward normal derivative which makes it an intrinsic 
definition, independent of the orientation chosen on the edge. 

2 

Definition 5.1. Cq' (Q) denotes the subspace of functions f in Co{G) which are twice 
continuously differentiable on Q°, and such that for every v ^ V and all I €z C the 
limit 

(5.2) fivi) = lim , no 

exists. Cq{Q) denotes the subspace of those functions f in Cq'^{Q) so that f" extends 
from Q° to a continuous function on Q. 

Thus Cl{g) consists of all / E Cq''^{Q) so that for every v £ V, the f"{vi), 
I G ^(^')^ are all equal. Assume that / G Cq{Q), and let v € V. The continuous 
extension of /" to v will simply be denoted by f"{v). Consider an edge I G C{v) 
incident with v. Then it is easy to see that f'{vi) exists (and is finite). On the other 
hand, in general for /, I' G C{v), I / we have f'{vi) / f'{vii). In other words, in 
general /' does not have a continuous extension from Q° to Q. Also, it is not hard to 
check that /' vanishes at infinity. 

The proof of the following lemma can be taken over with minor modifications from 
the standard literature, e.g., from [22, Chapter 6.1]. 



Lemma 5.2. For every Brownian motion X on the metric graph Q, the generator A of 
its semigroup U acting on Cq{Q) has a domain 'D{A) contained in Cq{Q). Moreover, 
for every f (^V{A), Af = 1/2 f". 

Consider data of the following form 

a = (a^, v£V)£ [0,1)^ 

(5.3) b={K„viGVc)G[0,l]^'^ 

c= (c„, veV)e [0,if 

subject to the condition 

(5.4) ay+ Ki + Cv — 1, for every v e V. 
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We define a subspace 'Ha,b,c of Cq{Q) as the space of those functions / in Cq{Q) 
which at every vertex v ^ V satisfy the Wentzell boundary condition 

(5.5) a.f{v) - Kf'M + \ cvf"{v) = 0. 

Now we can state and prove the analogue of Feller's theorem for metric graphs. 

Theorem 5.3. Suppose that X is a Brownian motion on a metric graph Q, and that 
T){A) is the domain of the generator A of its semigroup. Then there are a, b, c as 
in ^A\, so that V{A) = 'Ha,b,c- 



Remark 5.4. The case = 1, v £ V, would correspond to (zero) Dirichlet boundary 
conditions at the vertex v. The paths of the process associated with this boundary 
condition have to jump instantaneously to A when reaching the vertex, and by our 
requirement that the paths are right continuous this means that the process will never 



be at the vertex. But this is in contradiction to our assumption (cf. definition 3.11 
that the process with absorption at the vertex is equivalent to a Brownian motion with 
absorption in the endpoint (endpoints, resp.) of the corresponding interval. Therefore 



this stochastic process is not a Brownian motion on Q in the sense of definition 3.1 
and this case has to be excluded from our discussion. Also note that in this case the 
semigroup does not act strongly continuously on Co{Q), and therefore is in particular 
not Feller. 



The proof of theorem 5.3 has two rather distinct parts, and therefore we split it by 



proving the following two lemmas: 

Lemma 5.5. Suppose that X is a Brownian motion on a metric graph Q, and that 
!D(A) is the domain of the generator A of its semigroup. Then there are a, b, c as 
in ( |531 ), ^5A\ , so that C ?^a,fe,c- 

Lemma 5.6. Suppose that A is the generator of a Brownian motion X on Q with 
domain T>(A) C 'Ha,b,cfor some a, b, c as in ( |5.3| ), ( |5.4[ ). Then T>{A) = 'Ha,b,c 



Proof of lemma p3| Our proof follows the one in p2| Chapter 6.1] quite closely — 
actually, it is sufficient to consider a special case of the proof given there. 

We show that for every vertex v ^ V there are constants ay € [0, 1), G [0, 1], 
/ G C{v), Cv G [0, 1] satisfying ( |5.4| ), and such that all / in the domain T>{A) of the 



generator satisfy the boundary condition ( |5.5[ ). To this end, we let / G T^{A), fix a 



vertex v £ V, and compute Af{v). Let us consider the three cases for (3 mentioned 
in section |3] 

If /3 = 0, u is a trap, and Utf{v) = f{v) for all t > 0. Consequently, Af{v) = 0, 
and therefore 1/2 /"(f) = 0. Thus / satisfies the boundary condition ( |5.5| ) at v with 
= 0, c„ = 1, and = for all / G C{v). 

Next we consider the case where /3 G (0, +oo), i.e., v is exponentially holding. We 
know from the discussion in section[3]that then after expiration of the holding time the 
process jumps directly to the cemetery state. Therefore we get for t > 0, Utf{v) = 
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exp(— /3i)/(f ), and thus Af{v) + (3f{v) = 0, and the boundary condition ( |5.5| ) holds 
for the choice 



1 + 1 + /3' 



(5.6) = ., , ^ , = ^ ^ ^ , 6^,, = 0, Z G £(?;). 



Finally we consider the case that /3 = +c«, i.e., the process leaves v immediately, 
and in particular, v is not a trap. Therefore we may compute Af{v) in Dynkin's form, 
e.g., |[5} p. 140, ff.], fTS] p. 99]. As in ( |43| ) we let i?^,,, denote the hitting time of the 
complement of B^{v)7Then 



Recall the notation /;(e) for /(^) with ^ E ^ having local coordinates (/, e), / G i2, 
e E Then we get 

= Me)P,{X{H,,,) G + /(A)P,(X(i/,,,) = A) 

Ze£(»;) 

where the last equality follows from /(A) = 0. Let us denote 

= A) 



K{e) = l + r^ie) + € ^ n(e). 

The continuity of the paths of X up to the lifetime ( yields 

Y Pv{X{H,,,) G + Pv{X{H,,,) = A) = 1, 

and therefore equation ( |5.7| ) can be rewritten as 

lim[Af{v) + rA(.e)f{v)- ^ n(e) (//(e) - / W)) = 0. 

Since for all e > 0, K(e)~i < 1, it follows that 
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which by lemma [S!2| we may rewrite as 

Um(a.(6)/(^;) + ^c.(e)/»- ^ (e) ^'^'^ ~ ^^^^ ) = 0, 

where we have introduced the non-negative quantities 

rA(e) 



ay{e) 
c„(e) 



1 

en(e) 

K(6) 



Observe that for every e > 0, 

a^(e) + Ci,(e) + ^ 5^^(e) = 1. 

Therefore every sequence (e^, n G N) with e„ > and e„ | has a subsequence so 
that a„(e), c„(e) and Ki{e), I E ^(?^), converge along this subsequence to numbers 
o^, Cv, and h^^ respectively in [0, 1], and the relation (54 1 holds true. From the remark 
after definition l5.1l it follows that 

/Ke) - fjv) 

e 

converges with e | to /'(t;/), and therefore we obtain that for every vertex w G F, 
/ G T^{A) satisfies the boundary condition (5^1 with data a, b, c as in (53 1, (5^1. □ 

Before we can prove lemma [53| w e have to introduce some additional formalism. 

For given data a, b, c as in ( |5.3| ), ( |5.4| ), it will be convenient to consider 1-La,b,c 
equivalently as being the subspace of Cq'^(^) so that for its elements / at every 
V the boundary conditions ( |5.5| ) as well as the boundary condition 

(5.8) f"{vi) = f'ivk), for all I, k e C{v) 

hold true. Relation ( |5.8| ) is just another way to express that /" is continuous on Q. 

We consider the sets V, £, and Z as being ordered in some arbitrary way. With the 
convention that in C the elements of £ come first this induces also an order relation 
on C 

Suppose that / G C^''^{G). With the given ordering of £ and Z we define the 
following column vectors of length \£\ +2\Z\: 

fiV) = ((/e(0), ee£), {m, iez), {Mpi)ie z))\ 

nv) = [{I'M, ee£), im, i e ^))\ 

= ((/:'(0), eG£), (/r(0), i G X), ^ G X))*, 

where the superscript "f ' indicates transposition. 
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We want to write the boundary conditions ( |5.5| ), ( |5.8| ) in a compact way, and to this 
end we introduce the following order relation on Vc- For vi, v[, G Vc we set vi < v[, 
if and only if v ^ v' or v = v' and / < I' (where for V and C we use the order 
relations introduced above). For / as above set 

~f{V) = {f{vi), VI G Vc)\ 

fiV) = ifivi), VI G Vcf, 

fiV) = {f'ivi), VI e Vcf. 

Then there exists a permutation matrix P so that 

f{V) = PfiV), f'{V) = Pf'V), f"{V) = Pf"{V). 

In particular, P is an orthogonal matrix which has in every row and in every column 
exactly one entry equal to one while all other entries are zero. 

For every w G F we define the following |>C(v)| x matrices: 





/a„ 


•• 


• ^\ 











•• 


■ 






A{v) = 














1^0 


•• 


• oJ 










-K 




—bv, \ 

'\C(v)\ 


















B{v) = 














V 













/ 




/l/2c„ 








• 








1 


-1 





• 


■ 









1 


-1 


• 


• 




C{v) = 








1 


-1 • 


• 






I 








• 


• -V 





where we have labeled the elements in C{v) in such a way that in the above defined 
ordering we have h < h < ■ ■ ■ < ^C{v)\- Observe that C{v) is invertible if and only 
if Cv / 0. Define block matrices A, B, and C by 



A 



B 



C 



®C{v). 



Then we can write the boundary conditions ( |5.5[ ), ( |5.8[ ) simultaneously for all vertices 
as 



(5.9) 



Af{V) + Bf{V) + Cf"{V)=Q. 
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Consequently the boundary conditions can equivalently be written in the form 

(5.10) Af{V) + BfiV) + Cf'iV) = 0, 

with 



(5.11) 



A = P-^AP, B = P-^BP, C = P-^CP. 



We bring in the following two matrix-valued functions on the complex plane 

(5.12) Z±{k) = A±kB + K^C, KeC. 

Lemma 5.7. There exists R > so that for all k £ C with \k\ > R the matrices 
Z±{k) are invertible, and there are constants C, p > so that 

(5.13) \\Z±{k)-^\\ <C\k\p, \k\>R. 

Remark 5.8. The bound ( 5.13| ) is actually rather crude, but sufficient for our purposes. 



2(7) P 



Proof of lemma 5. 7 Since we have 

(5.14) Z±{k) = P-\A ± kB 

for an orthogonal matrix P, for the proof of the first statement it suffices to show that 
there exists R> such that 

A±kB + k^C 

are invertible for complex k outside of the open ball of radius R. For this in turn it 
suffices to show that for every vertex v the matrices 

A{v)±kB{v) + K^C{v) 








±Kbvi^ 








A 
















































-^2 


/ 



V 

are invertible for all «: G C with |k| > R. An elementary calculation gives 

,2 



dei[A{v) ± kB{v) + k^C{v)) = {a^±K ^ K + — 

iec(v) 



2\\C{v)\-l 



The choices k = ±1 together with condition ( |5.4| ) show that the polynomial of sec- 
ond order in k in the first factor on the right hand side does not vanish identically. 
Therefore, it is non-zero in the exterior of an open ball with some radius R^ > 0. 
Hence, we obtain the first statement for the choice R = max^gy R^. Moreover, from 
the calculation of the determinants above we also get for every v £ V and all k G C 
with \k\ > Ran estimate of the form 



(5.15) 



\det{A{v) ± kB{v) + K^C{v)) \ ^ < const. 
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Thus, using the co-factor formula for 

{A{v) ± kB{v) + k'^C{v))~^ 
we find with ( |5.15| ) the estimate 

\\{A{v) ± kB{v) + K'^C{v)y^\\ < C^\k\p% \k\ > R, 
for some constants Cy,py > 0. Consequently we get 

\\{A±KB + K'^Cy'^\\<C\K\P, \k\>R, 
for some constants C,p > 0, and by ( |5.14[ ) we have proved inequality ( |5.13| ). □ 
With these preparations we can enter the 



Proof of lemma pTS] Let the data a, b, c be given as in( |5.3| ), ( |5.4| ). We have to show 
that the inclusion 'D{A) C 'Ha,b,c is not strict. Let R = {Rx, A > 0) be the resolvent 
of A. Then for every A > 0, R\ is a bijection from Co{G) onto ^{A), that is, R^ is 
a bijection from 'D{A) onto Co{Q). 

Assume to the contrary that the inclusion T>{A) C 'Ha,b,c is strict. We will derive a 
contradiction. For A > consider the linear mapping Hx : f ^ Xf — 1/2/" from 
T~ia,b,c to Cq{Q). On 'C(A) this mapping coincides with , and R'^ is a bijection 
from ^{A) onto Cq{Q). Therefore our assumption entails that Hx cannot be injective. 
Hence for any A > there exists /(A) G 'Ha,b,c, /(A) / 0, with 

(5.16) Hxf{X) = Xf{X)-lf"{X) = 0. 

We will show that /(A) G 'Hafi.c satisfying ( |5.16| ) can only hold when /(A) = on 
G- It will be convenient to change the variable A to k = V2X, and there will be no 
danger of confusion that we shall simply write /(k) for /(A) from now on. Then the 
solution of ( |5.16| ) is necessarily of the form given by 

(5.17) fe{K,x) = re{K) 6 G X G M+, 

(5.18) /i(At,x) = r+(K)e''" + rr(K)e'^(^»-") i e I, x e [0, p^], 



and we want to show that for some k > 0, the boundary conditions ( |5.5| ) and ( |5.8| ) 
entail that re(K) = rf{K) = r~(«;) = for all e G i G T. For k > 0, define a 
column vector r{K) of length \£\ + 2|X| by 

r(K) = ((re(K), e G £), (r+(K), i G X), (r-(K), i G X))*, 
and introduce the (|<?| + 2|X|) x (|£^| + 2|X|) matrices 

X±{k) = 

— appropriately modified in case that f or X is the empty set — with the |X| x |X| 
diagonal matrices 

e'^'' = diag(e''''% iel). 
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Then the boundary conditions ( |5.5| ), ( |5.8| ) for /(k) read 

(5.19) Z{K)r{K) = 0, 
with 

(5.20) = {A + k'^C)X+{k) + kBX^{k). 

Thus, if we can show that for some k > the matrix Z{k) is invertible, the proof 
of the theorem is finished. Note that the matrix-valued function Z is entire in k, and 
therefore so is its determinant. Thus, if can show that k ^ det Z{k) does not vanish 
identically, then it can only vanish on a discrete subset of the complex plane, and for 
K in the complement of this set Z{k) is invertible. Write 

X±{k) = l±dX{K), 



with 

so that we can write 
with 



6X{k) 








Z{k) = Z+{k){1 + 6Z{k)), 



6Z{k) = Z+iK^ 6X{k). 
Observe that in case that X = 0, we obtain 6Z{k) = 0, and in this case the invertibility 



of Z{hc) for all hc with k > R follows from lemma 5.7 Hence we assume from now 
on that X / 0. Lemma [5?7] provides us with the bound 

< const. Ik]'', 

for all K & C,\k\ > R, and for some q > 0. On the other hand, we get 

\\6X{K)\\<e^P\ 

for all K < where po = minjgj pi. Therefore, there exists a constant i?' > so that 
for all K < —R' we have ||(5Z(k)|| < 1, and therefore for such k, Z{k) is invertible, 
i.e., det Z{k) ^ 0. Hence there also exists k > so that Z(^k) is invertible, and the 
proof is finished. □ 

Appendix A. Feller Semigroups and Resolvents 

In this appendix we give an account of the Feller property of semigroups and re- 
solvents. The material here seems to be quite well-known, and our presentation of 
it owes very much to [32|, most notably the inversion formula for the Laplace trans- 



form, equation ( |A.3| ) in connection with lemma A.6 On the other hand, we were 
not able to locate a reference where the results are collected and stated in the form 
in which we employ them in the present paper. Therefore we also provide proofs for 
some of the statements. 

Assume that (E^ d) is a locally compact separable metric space with Borel a- 
algebra denoted by B{E). B{E) denotes the space of bounded measurable real valued 
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functions on E, Cq{E) the subspace of continuous functions vanishing at infinity. 
B{E) and Cq{E) are equipped with the sup-norm || • ||. 
The following definition is as in [33] : 

Definition A.l. A Feller semigroup is a family U = {Ut, t > 0) of positive linear 
operators on Cq{E) such that 

(i) = id and \\Ut\\ < 1 for every t > 0; 

(ii) Ut+s = Ut° Us for every pair s,t> 0; 

(iii) lim^o \\Utf - /II = for every f G Co{E). 

Analogously we define 

Definition A.2. A Feller resolvent is a family R = {R\, A > 0) of positive linear 
operators on Cq{E) such that 

(i) II^aII < for every A > 0; 

(ii) R\ — Rfj, = {fi — X)Rx o Rfifor every pair \ fi > 0; 

(iii) liniA^oo ||Ai?A/ - /II = for every f G Co{E). 

In the sequel we shall focus our attention on semigroups U and resolvents R asso- 
ciated with an E'-valued Markov process, and which are a priori defined on B{E). 
(In our notation, we shall not distinguish between U and R as defined on B{E) and 
their restrictions to Cq{E).) 

Let X = {Xt, t > 0) be a Markov process with state space E, and let {P^, x G 
E) denote the associated family of probability measures on some measurable space 
{Q.,A) so that Px{Xq = x) = 1. Ex{-) denotes the expectation with respect to P^. 
We assume throughout that for every / G B{E) the mapping 

{t,x)^Ex{f{Xt)) 

is measurable from M+ x E into M. The semigroup U and resolvent R associated 
with X act on B{E) as follows. For / G S(^), x G -E, i > 0, and A > set 

(A.l) Utf{x) = Ex{f{Xt)), 

/•oo 

(A.2) Rxf{x)= / e~^'Utf{x)dt. 

Jo 



Property (i) of Definitions|A. l|and|A.2 is obviously satisfied. The semigroup property. 



(ii) in Definition A. 1 follows from the Markov property of X, and this in turn implies 



the resolvent equation, (ii) of Definition A.2 Moreover, it follows also from the 



Markov property of X that the semigroup and the resolvent commute. On the other 
hand, in general neither the property that U or R map Co{E) into itself, nor the strong 



continuity property (iii) in Definitions A. 1 A.2 hold true on B{E) or on Cq{E). 



If 1^ is a subspace of B{E) the resolvent equation shows that the image of W 
under Rx is independent of the choice of A > 0, and in the sequel we shall denote 
the image by RW. Furthermore, for simplicity we shall write UCq{E) C Cq{E), if 
Utf G Cq{E) for alH > 0, / G Cq{E). 

Theorem A.3. The following statements are equivalent: 
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(i) U is Feller. 

(ii) R is Feller. 

(iii) UCoiE) C CoiE), and for all f G Co(^), x e E, lim^o Utfix) = f{x). 

(iv) UCq{E) C Co{E), and for all f G Co{E), x e E, limx^oo ^Rxf{x) = 
fix). 

(v) RCo{E) C Co{E), and for all f G Co{E), x £ E, lim^o Utf{x) = f{x). 

(vi) RCo{E) C CoiE), and for all f G Co(^), x G E, liuix^oo >^R\fix) = 
fix). 

We prepare a sequence of lemmas. The first one follows directly from the domi- 
nated convergence theorem: 

Lemma A.4. Assume that for f G B{E), Utf ^ f as t I 0. Then \R\f f as 
A — > +00. 

Lemma A.5. The semigroup U is strongly continuous on RB[E). 

Proof. If strong continuity at t = has been shown, strong continuity at t > follows 
from the semigroup property of U, and the fact that U and R commute. Therefore it 
is enough to show strong continuity at t = 0. 

Let / G B{E), A > 0, t > 0, and consider for x G -E the following computation 

UtRxfix) - Rxfix) 

roo /'OO 

= / e-^'E4f{Xt+,))ds- e-^'E,{f{Xs))ds 
Jo Jo 

poo poo 

= e^* / e'^'E, {f{Xs)) ds - / e'^'E, (/(X,)) ds 
Jt Jo 

poo pt 

= (e^* - 1) / e-'^'E^ifiXs)) ds - / e-""' E^{f{Xs)) ds 
Jt Jo 

where we used Fubini's theorem and the Markov property of X. Thus we get the 
following estimation 

poo pt 
^Xt '\ \ I As I / As 



\\UtRxf-Rxf\\ < (^(e-^ -1) e-^'ds + e-^' ds 
= f 11/11, 

which converges to zero as t decreases to zero. □ 
For A > 0, t > 0, / G B{E), x £ E set 

(A.3) U.^fix) = V ^ \ n\ e"^* Rnxf{x). 

n=l 

Observe that, because of nA||i?„A/|| ^ ||/||, the last sum converges in B{E). 
For the proof of the next lemma we refer the reader to |32[ p. 477 f]: 
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Lemma A.6. For all t > 0, / G RB{E), U^f converges in B[E) to Utf as A tends 
to infinity. 

Lemma A.7. IfUtCo{E) C Co{E)forallt > 0, then RxCq{E) C Cq{E), for all 
A > 0. If RxCq{E) C Cq{E), for some A > 0, and RxCq{E) is dense in Co{E), 
then UtCo{E) C Cq{E) for all t > 0. 

Proof Assume that UtCoiE) C Cq{E) for all t > 0, let / G Co(£^), x G and 
suppose that n G N) is a sequence converging in {E, d) to x. Then a straightfor- 
ward application of the dominated convergence theorem shows that for every A > 0, 
R\f{xn) converges to Rxf{x). Hence Rxf G Cq{E). 

Now assume that that RxCq{E) C Cq{E), for some and therefore for all A > 
0, and that RxCo{E) is dense in CoiE). Consider / G RCq{E), t > 0, and for 
A > define U^f as in equation ( |A.3| ). Because Rnxf G Co{E) and the seri es in 



formula ( A.3 1 converges uniformly in x G E, we get U^f G Cq{E). By lemma |A. 6 
we find that U^f converges uniformly to Utf as A — +oo. Hence Utf G Cq{E). 
Since RCq{E) is dense in Co{E), Ut is a contraction and Cq{E) is closed, we get 
that UtCo{E) C Co(^) for every f > 0. □ 

The following lemma is proved as a part of Theorem 17.4 in [20 1 (cf. also the proof 
of Proposition 2.4 in |[33|). 

Lemma A.8. Assume that RCq{E) C Cq{E), and that for all x e E, f e Cq{E), 
limA-j>oo ^R\f{x) = f{x). Then RCq{E) is dense in Cq{E). 

If for all / G Co{E), x ^ E, Utf{x) converges to f{x) as t decreases to zero, 



then similarly as in the proof of lemma |A!4 we get that XRxf{x) converges to f{x) 



as A — )• +00. Thus we obtain the following 

Corollary A.9. Assume that RCq{E) C Cq{E), and that for all x £ E, f £ Cq{E), 
lim^o Utf{x) = f{x). Then RCq{E) is dense in Cq{E). 

Now we can come to the 



Proof of theorem A.3 We show first the equivalence of statements (i), (ii), (iv), and 
(vi): 



"(i)=^ (ii)" Assume that U is Feller. From lemma A.7 it follows that RxCq{E) C 
Co{E), A > 0. Let / G Cq{E). Since U is strongly continuous on Cq{E), lemma A.4 
implies that XRxf converges to / as A tends to +oo. Hence R is Feller. 

"(ii)=^ (vi)" This is trivial. 



"(vi)=^ (iv)" By lemma A.8[ RCq{E) is dense in Co{E), and therefore lemma A.7 



entails that UCq{E) C Cq{E). 



"(iv)=^ (i)" By lemmas |A77 and |A.8 RCq{E) is dense in Cq[E), and therefore by 



lemma A.5 U is strongly continuous on Cq{E). Thus U is Feller. 

Now we prove the equivalence of (i), (iii), and (v): 
"(i)=^ (iii)" This is trivial. 
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"(iii)=^ (v)" This follows directly from Lemma A.7[ 

"(v)=^ (i)" By corollary A.9[ RCq{E) is dense in Cq{E), hence it follows from lem- 



ma 



A.7 that U Co{E) C Cq{E). Furthermore, lemma A.5 implies the strong continu- 

□ 



ity of U on RCq{E), and by density therefore on Cq{E). (i) follows. 
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